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We consider a dc-driven damped sine-Gordon model with a small nonlinear spatial-disorder term, 
onto which a sinusoidal modulation is superimposed. It describes, e.g., a weakly disordered system 
with a regular grain structure. We demonstrate that, at the second order of the perturbation 
theory (with respect to the weak spatial disorder), the periodically modulated disorder gives rise 
to an effective periodic potential. Dynamics of a kink moving in this potential is studied in the 
overdamped limit, using the adiabatic approximation, the main objective being to consider depinning 
of a trapped kink. The analytical results are compared with direct dynamical simulations of the 
underlying model, as well as with numerical results using the collective-coordinate approach but 
without the mean-field approximation. It is found that a critical force for the depinning of a kink 
trapped by the periodically modulated weak spatial disorder is much larger than that predicted by 
the mean-field approximation. 

PACS: 05.45.Y, 61.43.Bn, 71.55.Jv 
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I. INTRODUCTION 

The interplay of nonlinearity and disorder is a well- 
known fPRW in theoretical and experimental studies of 
solitonsETcfB. A paradigm one-dimensional model which 
has numerous physical applications and, in particular, 
allows one to consider interaction of topological solitons 
(kinks) with a random spatial inhomogeneity is a per- 
turbed sine-Gordon (SG) equation. In-,the simplest case, 
the corresponding model takes a formElu, 

- ®xx + [1 + e(x)\ sin$ = -(r/8)$ t - /, (1) 

where $ is a dynamical order parameter, T is a dissipa- 
tion constant, / = const represents a driving force (exter- 
nal dc field), and the disorder function e(x) is determined 
by its white-noise Gaussian correlations, 



(e(x))=0, (e(x)e(x )) = V 2 6(x - x ), 



(2) 



where r\ is the disorder amplitude. In this work, the 
perturbation coefficients rj 2 and / will be treated as small 
parameters, while Y is not necessarily small. Note, that 
similar disordered models can be defined not Qaly in the 
continuum form, but also on aperiodic latticesO. 

In some systems, the spatial disorder can be subject 
to a superimposed regular (periodic) modulation. For 
instance, this may be the case if the underlying slightly 
disordered physical medium has a regular grain structure. 
Then, the model ([!]) is modified as follows, in the simplest 
case when the superimposed modulation is sinusoidal, 



*« - + [1 + e(x) sm(kx)} sin $ = -(r/8)$ t - /, 



(3) 



2ir/k being the modulation period. 

As concerns the form of the model (||), it is relevant 
to mention that, for any finite k, one can replace the 
modulation function sin(kx) by cos(fcir), shifting x — > 
x + 7r/2fc. If the system is very long and e(x) obeys the 
Gaussian correlations (||), the modulation functions are 
tantamount to each other, except the case k — (or 
\k\<l/l, if the system has a finite length 1). Actually, 
for k — one comes back to the model p) with the 
unmodulated spatial disorder. 

Dynamics of topological nonlinear excitations in the 
form of kinks in the model (||) is a subject of the present 
work. In the lowest-order approximation, the random 
perturbation does not induce any effective potential for 
the kink. However, we demonstrate in section 2 that 
an effective periodic potential is generated at the second 
order of the perturbation theory. It is interesting that the 
amplitude of the effective potential vanishes at k 2 = 1/2. 

In section 3 we study motion of the driven kink in 
this effective potential. The corresponding equation of 
motion is derived in the adiabatic approximation (treat- 
ing the kink as a quasiparticle). This approximation is a 
version of the well-known collective-coordinate technique, 
which was applied to the (weaklyl-disordered models in 
various contexts (see, e.g., Refs.BQ), allowing one not 
only to derive an effective equation of motion, but also ' 
study dynamical statistical characteristics of the kin! 
We focus on the case when the overdamped approxima- 
tion applies to the description of the kink motion (i.e., 
the inertia term in the corresponding equation of motion 
is neglected). In its direct form, the latter approximation 
assumes that the dissipation is stronger than the disor- 
der, or, in terms of Eq. ([!]) r ^> 77, in the generic case 
with fc ~ 1. However, in this work we actually focus on a 
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transition between the pinning and free-motion dynam- 
ical regimes (i.e., depinning of the trapped kink) with 
the increase of the driving force /. It is quite reasonable 
to assume that the inertia plays a negligible role in the 
depinning process, hence the overdamped approximation 
can be applied in a broad parametric region. Note that 
depinning of kinks trapped by the unmodulated spatial 
disorder in the model (|lj) was earlier studied in detail in 
RefM 

Section 4 displays results of numerical simulations of 
the model, which consists of two parts: direct dynami- 
cal simulations, and numerical analysis of the collective- 
coordinate description without using the mean-field ap- 
proximation. The latter approach yields accurate results 
for the critical (minimum) force / cr necessary to depin 
the trapped kink. The critical force turns out to be much 
larger than that predicted by the mean-field approxima- 
tion. In particular, / cr is predicted by the mean-field 
approximation to scale ~ rj 2 (see Eq. (g)), while the ac- 
tually found critical force scales ~ 77, the dependence on 
the modulation wave number k (see Eq. also being 
quite different. 



II. SINE-GORDON KINK IN A PERIODICALLY 
MODULATED RANDOM POTENTIAL 

Equation (J3J) can be solved by iterations, assuming a 
small amplitude of the noise. The unperturbed kink so- 
lution is 

$o(x - £) = 4tan _1 exp(x - £) , (4) 

£ being the coordinate of the kink's center. 

A first-order correction to the kink waveform, <j>\ = 
<&(af) — $0j satisfies the equation 

(0i)tt - (01 ) axe + cos $o0i = -e(x)sin(fcx)sin$ . (5) 

To solve Eq. (g), we need its Green's function, which is 
defined in a standard way by means of the equation, 

{§2 ~ -H~2 + V"[*o(x)] } G(x, x', t) = 8{x - x')5(r), 

(6) 

where r = t — t' and V($>) = 1 — cos$ is the SG potential. 
To simplify the formulas, we temporarily set £ = 0. 

The necessary Green's function has-been found in an 
explicit form by Flesch and TrullingerQ (see Eq. (3.29)): 

G sg {x,x',t) 



[6{t - \z\)/2] {J (s) + 2 [sJi{s) - 2tA 1 {w, s)] 

-/3 3 sgn(z)[-(r + \z\)J (s) + 2rA ( W , s)}} . (7) 



Here, z = x — x 1 , s = yr^ 



t — \z\, 6(t) is the 



1, /?3 = tanh(x') — tanh(x). Further, Jq(s) and Ji(s) are 
the Bessel functions and Aq(w,s) and Ai(w,s) are the 
modified Lommel functions defined as 



, / w \2m+n 
A n (w,s) = y — J J2m+n{s) 



(8) 



m=0 




FIG. 1. The effective potential U eS {k,£). 
Then, <f>\ can be obtained in the form 



4>i(k,x) 



dx C(x') sin(<po(x' — £)) 
x sm(kx')G(x — x') , 



(9) 



where we have restored the kink's coordinate £. 

The perturbation on the right-hand side of Eq. (||) gives 
rise, in a straightforward wayo, to an effective potential 
for the kink treated as a quasiparticle: 



U e s(k,£)= / dx e(x) sin(kx)4>i(k, x) sin <&o{x — £). 



(10) 

Then, inserting the expression (^) for 4>\{k, x) into Eq. 
(|To|), one obtains 



00 00 



dx / dx e(x)e(x ) sin $0(2; — 



-00 —00 



xsin$o(x — sm(kx) sin(kx )G(x — x ). (11) 



Heaviside's step function, and 02 = tanh(a;) tanh(x') 



The potential (|llj) can be averaged with regard to 
(||). Using the elementary expressions for the SG kink, 
cos $0(2;) = 1 — 2sech 2 x and sin<f>o = 2 sinhxsech 2 a;, 
yields the average potential 
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FIG. 2. The phase diagram for r\ = 1. The regions of the 
free-motion (|A/_B| < 1) and pinning > 1) regimes 

are, respectively, white and grey. At the boundary between 
them, \A/B\ = 1. 
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cosh 4 [x 

7rfc(l 



sin 2 (fez) 



2k 





2 ^ cos(2fc£) 
sinh(7rfc) 



(12) 



Evidently, the average potential (U e g(k, £)) ( |12[ ) is pe- 
riodic. It is plotted in Fig. 1 as a function of £ and 
k. 

Minima of the potential (|l2| ) are at the points £ m i n = 
(<n-/2k) (2n + 1), n = 0,1,..., for fc 2 < 1/2, and 
^min = (n/k)n for A; 2 > 1/2. Close to the minima, 
£ = Cmin + <5, kS « 1, the potential gives rise to 
small nearly harmonic oscillations with the frequency 
n 2 = 477 2 7rfc 2 |2fc 2 - 1\/ (3sinh7rfc). 



III. DRIVEN MOTION OF THE KINK IN THE 
EFFECTIVE POTENTIAL: ADIABATIC 
APPROXIMATION 

In order to investigate uniformly driven damped mo- 
tion of the kink in the effective periodically modulated 
random potential (|l2|), we derive, in the framework of 
the standard adiabatic (quasi-particle) approximation^, 
an equation of motion for the kink's center, 



dt 2 



dt 



2irf- 



4n 2 irk 2 (l - 2fc 2 ) 
3 sinh(7rfc) 



sin(2fc£) (13) 



(here, the "nonrelativistic" approximation, (d^/dt) <C 1, 
is additionally adopted; nate that the kink's mass in the 
present notation is m = 80). 

As it was explained above, we will chiefly be interested 
in the depinning threshold, when the inertia term in Eq. 
( |l3| ) may be neglected. Then, the equation of motion 
simplifies to 



d£, _ 2tt/ 4?y 2 7rfc 2 (l - 2fc 2 
~dl ~ ~V~ 3r sinh(Trfc) 



■ sin(2fc£) 



B + Asin(2fc£), 
(14) 



where A = 4r] 2 irk 2 (l - 2fcf)/ [3r sinh(Trfc)] and B = 
2irf/T. Comparing Eqs. ( ^3| ) and (0), one can verify 
that the inertia term is indeed negligible provided that 
r ^> \A\k. In the general case (k ~ 1), this simply means 
r> i), i.e., the dissipative constant must be essentially 
larger than the amplitude of the spatial disorder. 

A solution to Eq. (jlj) is obvious. One has to dis- 
tinguish the cases (i) \A\ = \B\, (ii) \A\ < \B\, and (iii) 
\A\ > \B\. The driving and potential terms dominate, re- 
spectively, in the regions (ii) and (iii) , while (i) represents 
a border between them, characterized by the value 



/ b = (2/fcV/3)|l-2fc 2 |/sinh(7rfc) 



(15) 



of the driving field. The corresponding phase diagram is 
plotted, in terms of the underlying parameters k and /, 
in Fig. 2, where the boundary between the regions (iii) 
(grey) and (ii) (white) corresponds to the critical field 
( |l5| ) . Basic features of these regimes, that can be deduced 
in an analytical form from Eq. (|14|), are summarized 
below. 
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FIG. 3. Free (progressive) motion of the kink, shown in 
terms of x(t) (a) and velocity v(t) (b). Here and in the figures 
below the results are displayed for F = 10 and V — 5 by the 
thin and thick lines. The other parameters are f] = 1, k = 1, 
and / = 0.1. 

(i) The border regime, \A/B\ = 1: Eq. ( |l4| ) gives 
rise to semi-stable fixed (stationary) points, £ s = {An — 
l)(-7r/4fc), n = 0, ±1,..., which seem like attractors and 
repellers to the left and right of them, respectively. In 
this case, an analytical solution to Eq. ( |l4| ) takes an 
implicit form, 



cos(2fc£) 



2kA{l + sin(2fc£)) 



= t-t c 



(16) 



to being an initial moment. 

(ii) The free (progressive) motion regime, < 1. 

In this case, Eq. ( |l4| ) does not have stationary solutions. 
A time-dependent solution is 
tan (fc£(i)) 



-A/B+ y/\ -A 2 IB 2 tan k\J B 2 — A 2 (t - t ) 



(17) 

A typical example of this solution is displayed in Fig. 
3. It represents a superposition of a linear trajectory and 
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FIG. 4. The average velocity v vs. the driving field for 
k = 1, T) = 1. The phase transition occurs at ft = 0.0577264. 
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FIG. 5. The average velocity v vs. the modulation wave 
number k for / = 0.04, r\ = 1. 



periodic oscillations. The linear trajectory corresponds 
to the average velocity v = lim h [£(i + T) — £(to)] = 



V-B 2 — A 2 . In a neighborhood of the free-motion's bor- 
der, at / — fb — > +0, one has 



v » (27T/T) V2/b(/-/b). 



(18) 



The average- velocity dependence as per Eq. (|18|) exhibits 
a typical mean-field critical (cusp) behavior close to the 
border, v playing the role of an order parameter (Figs. 4 
and 5). 

The time-periodic component of the velocity (see Fig. 
3b) is characterized by the frequency Li — 2fcV B 2 — A 2 , 
whose critical behavior at /— /& —> +0 is the same as that 
of the velocity ([l8]). The mean-square deviation of the 

T 

trajectory from the linear one, A = lim ^ J dt(^(t) — 

T^oo o 

_ _ T 

vt—Q 2 , where C = lim A J dt(£(t)—vt), exhibits too the 

T^oo 

critical behavior as function of /. The /-dependencies of 
the reduced frequency uj/2tt and of the deviation A are 
shown in Fig. 6. 

Note that, representing the kink's law of motion in the 
form £ = vt + C(~t)i where is a time-periodic function 
with a small amplitude, one can expand the kink solu- 
tion as follows (here, the "nonrelativistic" approximation 



<C 1 is assumed): 

<f)(x,t) = 4tan _1 {exp[a:- (vt + ((t))}} 
4 tan _1 {exp(a; — vt)} + 2((t)sech(x — vt). 



(19) 



The small oscillating term in Eq. (|l9|) may be a 
source for emission of radiation, provided that it con- 
tains spectral components with frequencies u> > 10. The 
Fourier decomposition of the source, defined as £(t) = 

CXD 

Y^, a n exp(inwi), is illustrated by Fig. 7. Evidently, 

n=— oo 

the density of the spectral components decreases with 
the distance from the boundary between the pinning and 
free-motion regions, i.e., with the increase of /. It is also 
obvious that the intensity of the generated radiation is 
extremely small, as the maximum frequency shown in 



Fig. 7 is « 0.377 of the threshold value, ui = 1, beyond 
which the emission of radiation takes place. 

(hi) The pinning regime, > 1. In Fig. 5, one 

has the pinning (v — 0) and free-motion (v ^ 0) regions 
of k in accord with the diagram in Fig. 2. In the pinning 
regime, Eq. (|l4|) exhibits two sets of fixed points: po- 
tential minima, £ a = {%/k)n — (l/2fc) sin -1 (B/A) , n = 
0, ±1, which are attractors (stable fixed points), and 
potential maxima, which are repellers (unstable points), 
C r = (4rc - l)7r/4fc + (l/2fc) sin" 1 (B/A). In the pinning 
regime, a time-dependent solution to Eq. (jl4|) is (cf. Eq. 

tan (fcf (*)) 



-A/B+ ^A 2 /B 2 - ltan ky/ A 2 - B 2 (t - t ) 



(20) 

This solution describes overdamped motion of a parti- 
cle trapped in a well of the tilted potential (U e g). The 
stationary state is reached after a short transient. 



IV. NUMERICAL SIMULATIONS 

A. Direct simulations of the perturbed sine-Gordon 
equation 

To get direct insight into the dynamics of the present 
weakly disordered model and check the above analytical 
results, we solved Eq. (0) numerically. For this purpose, 
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FIG. 6. Frequency of the oscillating part of the velocity and 
the mean square deviation of the trajectory from the linear 
one in the free-motion regime (ii), u(f) (a) and A(/) (b) for 
r] = 1, k = 1. 
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FIG. 7. Fourier spectra of the periodic part of the tra- 
jectory for r) = 1, k = 1, and V = 10: / = 0.057942 (a), 
/ = 0.076369 (b), and / = 0.115665 (c). 

Eq. (|h was transformed to a set of difference equa- 
tions in the spatial coordinate (differential in the time 
coordinate) and solved, by using the eighth-order explicit 
Runge-Kutta schemeEj, with a stepsize control such that 
the time step was dynamically changed within the range 
from 0.05 to 0.3. It has been chosen so that the relative 
(per one site) error at each step did not exceeded 10~ 10 . 
In order to study the kink dynamics within rather big 
time intervals (for our computations the required time 
lap was < t < 10 4 ) we have used a shifting computa- 
tional domain with boundaries which absorbed the out- 
going waves. 

At t = 0, we took, for the initial conditions, the well- 
known approximate (perturbative) static Junk solution 
valid in the absence of the spatial disorders, 

$k(s) =<p s +4tan- 1 (exps), s= - Z^= , (21) 



Uo 2 = (i + (r/27r/) 2 )- 



where 



sin" 1 / 



(22) 



takes into account the shift of the background under the 
action of the dc drive. 

In the simulations we used 1000 different realizations 
of the disorder {e} and then averaged the results over the 
realizations with equal weights. 

In Fig. 8, we compare the kink velocity, as obtained 
from the direct simulations of Eq. (||), and those pre- 
dicted by Eq. (|l4l) . The mean field velocity overestimates 
the exact one what is a typical mean field feature. 

The numerically found time dependence of the kink's 
velocity in the localization regime is plotted in Figs. 9 
and 10. The corresponding dependence of the inverse 
total distance travelled by the kink, in this regime, during 
the simulation time is shown, as a function of the driving 
force, in Fig. 11. It is evident from this figure that, 
while with small / the kink quickly gets trapped, with 
the increasing / it keeps moving for a much longer time. 



FIG. 8. Comparison of the velocity-vs.-time curves ob- 
tained from the numerically exact (diamonds) and adiabatic 
(stars) (Eq. (0)) solutions for / = 0.09, t) = 0.1, k = 1, and 
T = 2. 

When approaching the critical point, the trapping time 
enormously increases, therefore the direct numerical de- 
termination of the critical force / cr from the dynamical 
simulations is difficult. Besides, the calculations are con- 
strained to a finite size of the system, hence the critical 
field would be size dependent. Moreover, at low T, when 
//r > 1, the one-kink model fails, as nucleation of kink- 
antikink pairs sets in with the accumulation of energy. 
Therefore, investigation of the critical region (transition 
between the trapped and freely moving kinks is not really 
possible on the basis of the direct simulations from Eq. 



B. Numerical simulations of the collective 
coordinate equations 

Assuming low fields /, we will now follow the Rice s 
collective coordinate approacht 2 ], based on a modified 
ansatz for the kink, cf. Eqs. (f|) and ( [H| ) 



(x, t) = (f> s + 4 tan 1 exp 



L(t) 



(23) 



The difference from Eq. (|J) is that the present ansatz 
admits a variable kink's width L(t). Inserting (^3|) into 
the model's Hamiltonian 

oo 

H= J dx{(<$> 2 t +<P 2 x ) /2+[l + e{x)sin{kx)} 

— OO 

x(cos0 s - cos$) - fx$ x }, (24) 

one arrives at the effective Hamiltonian! 
L 3L 

H cS = -p\ + -^P 2 L + U(L) + V({e}, L, $ + 2irfZ 

(25) 

for the collective coordinates £, L and the canonically 
conjugate momenta 
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FIG. 9. Evolution of the soliton velocity for different values 
of the driving / = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7. V = 2, k = 1, 
r\ — 0.1 and initial position xo = 12.5. 



8 d£ 



2tt 2 dL 



Pi Ldt' PL 3L dt 



(26) 

In the effective Hamiltonian (|25|), 

U(L) = 4L- 1 + 4cos(0 s )L (27) 
is the effective potential in the case without disorder, and 

oo 

V({e},L,£) = / dxe(x)sm(kx){cos<j) s - cos$) (28) 



is an effective random potential. 

Taking into account the friction, we obtain the fol- 
lowing system of equations of motion for the two kink's 
degrees of freedom: 



dpL 
dt 



3p7 



4ir 2 



P\ d_ 
16 + dL 



Figs. 12 and 13 display comparison of the kink's veloc- 
ities and widths, as obtained from the direct simulations 
of Eq. (|3|) and from simulations of the collective coordi- 
nate equations (^9|) and (^fj|). A satisfactory agreement 
between the results of the two approaches is evident. 

Eq. (^9j) makes it possible to determine the critical force 
/ cr necessary to detrap the kink. For a single potential 
well (trap), / cr is simply equal to the maximum value of 
the local force induced by the trapping potential. How- 
ever, a spatially random potential assumes the presence 
of a collection of the wells with randomly distributed pa- 
rameters. Strictly speaking, one cannot uniquely define 
a global value of / cr in this case. Instead, it is necessary 
to consider a dynamical problem of gradually datrapping 
the kinks with the increase of the driving forceS. 
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FIG. 10. The same as in Fig. 9, for r) = 0.01, 0.1, 
/ = 0.04, r = 2, k = 1, and initial position xo = 12.5 for 
5.10 3 realizations of the random potential. 

Nevertheless, it is possible to give a reasonable defi- 
nition of the critical detrapping force as that which is 
equal to the local maximum value of the spatially ran- 
dom force, averaged over the length of the system (i.e., 
over all the local maxima): 



(31) 



2nf CI ={-V({e},0 



the angular brackets standing for averaging. Here, the 
kink's width in the potential V (^8j) should be set equal 
to the unperturbed value L = 1. 

If we suppose equal statistical weights for various local 
extrema of the-ijunction dV/ dL. then the averaging theo- 
rem from Ref.Eil (used e.g. inQ) can be applied directly, 
yielding 



(29) 



3jV( W ,*)) 



7T M 2 



Ma 



(32) 



[U(L) + V({e},L,C)]=0. (30) 




FIG. 11. Dependence of the inverse mean trajectory A 1 
on the force / for r = 2, fe = 1, rj = 0.1 and initial position 
x = 12.5. 
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20 40 60 80 100 

FIG. 12. Comparison of the averaged velocities vs. time 
from 

the direct numerical simulations of Eq. (^) (squares) with 
that from the numerical solution of the collective coordinate 
equations (diamonds) for tj = 0.1, / = 0.08, k = 1, T = 2. 



where, in the present case, 
' 'd 2 V 



Mi = 



Ma 



= V / dxshf(kx)q (z), (33) 



2 v oc 

(if) ) 1,/J f J*™ 2 «>:r)q (:). CUI 



and 2q , (z) = (cos c/f> s ) sech 2 (z) — (sin^,) sinh(z)sech 2 (,2:) 
(with z=(x- Q/L). 

The integrals in ( p3|) and (|^) can be evaluated in a 
straightforward way, 



Ij = (l-fl)Q +flR 3 , 3 = 1,2; 



(35) 



where 



32 / o 2 1 

h = 2i - 16cos(2fc05(7rfcL) f -Pi - -P 2 + — P 3 



(36) 



Pi = | - COB(2fc05(7Tfci) (l - |p + f P 2 - |P 3 



(37) 



Q 2 = ^ - 128cos(2^)S(7rfcL) Qpi - ^P 2 



^T7 p 4 , (38) 



4 2 
"2T 3 _ 315 Pi 



P 2 = ^-co S (2 W (^P)fl-fp 1 + fp 2 
15 \ 3 15 

256 „ 2.24 2 n \ /on . 

-^r P3 + ^35 P4 )■ (39) 
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FIG. 13. The same for the averaged widths at the same 
values of parameters as in Fig. 13. 

n 

where S(x) = x/sinhx ,P„ = II (( kL ) 2 + l2 )- Th en, 

i=i 

for /cr, one arrives at expression 



fcr 



n 



(40) 



The resulting function f CI (kL,£) is plotted in Fig. 14. 
To compare it with the mean-field result shown above 
in Fig. 2, we have set here £ = 0. A drastic difference 
between the two results is evident. Note also that the 
linear dependence of / cr on the disorder strength r/ is in 
contrast to the quadratic dependence on r\ in the mean 
field result (§J). 

It is relevant to stress that a characteristic critical de- 
pinning force, which was founcLfor the model ([!]) with the 
unmodulated disorder in Ref.S, also scaled linearly with 
1}. In other words, the depinning process is dominated 
by the underlying spatial disorder, rather than by the 
fact that the disorder is subject to a regular modulation. 
Thus, the mean-field approximation (alias, straightfor- 
ward perturbation theory for the weak disorder) can not 
adequately describe depinning. Nevertheless, the modu- 
lation gives rise to a dependence of the pinning threshold 
upon the modulation period, see Fig. 14. In the limit of 
large k, the oscillations are dominated by the random- 
ness and the case is reduced to that of the pure random 
noise. However, the most interesting is the region of the 
length scale competition hr 1 ~ L < 1. 

Note, that a strong nonperturbative response of soli- 
tons to weak random fluctuations in a narrow vicinity 
of a transition between different states is also known-, in 
other models, e.g., a dual-core nonlinear optical fiberu. 



V. CONCLUSION 

In this work, we have introduced a dc-driven damped 
sine-Gordon model with a small nonlinear spatial- 
disorder term, on which a sinusoidal modulation is super- 
imposed. It describes, e.g., a weakly disordered system 
with a regular grain structure. It was demonstrated that, 
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in the mean-field approximation (second order of the per- 
turbation theory), the periodically modulated disorder 
gives rise to an effective periodic potential. Dynamics of 
a kink moving in this potential was studied in the over- 
damped limit, using the adiabatic approximation, the 
main objective being to consider depinning of a trapped 
kink. The analytical results were compared against direct 
numerical simulations of the underlying model, as well as 
against numerical results using the collective-coordinate 
approach but without the mean-field approximation. It 
has been found that depinning of a kink trapped by the 
weak spatial disorder is, in reality, drastically different 
from what predicted by the mean-field approximation. 



tcr (K) 




FIG. 14. Critical field / cr vs. wave number of the modu- 
lating field k. Full line is for £ = 7r/2fc, dashed for £ = 0. 
Quantitative as well as the qualitative difference to the mean 
field version in Fig. 2 is evident. 
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